arXiv:1506.08019vl [math.OC] 26Jun2015 


MULTIOBJECTIVE APPROACH TO OPTIMAL CONTROL 
FOR A DENGUE TRANSMISSION MODEL 


ROMAN DENYSIURi 
roman.denysiukOalgoritmi.uminho.pt 

HELENA SOFIA RODRIGUES^’^ 
sofiarodriguesQesce.ipvc.pt 

M. TERESA T. MONTEIROi’^ 
tmOdps.uminho.pt 

LINO COSTAI’4 

lacOdps.uminho.pt 

ISABEL ESPIRITO SANTO^’^ 
iapinhoOdps .uminho.pt 

DELFIM F. M. TORRES^’^ 
delfim@ua.pt 


^Algoritmi R^D Center, University of Minho, 4710-057 Braga, Portugal 
^School of Business Studies, Viana do Castelo Polytechnic Institute, 4930-678 Valenga, Portugal 
^CIDMA R&D Center, University of Aveiro, 3810-193 Aveiro, Portugal 
“^Department of Production and Systems Engineering, University of Minho, 4710-057 Braga, Portugal 
^Department of Mathematics, University of Aveiro, 3810-193 Aveiro, Portugal 


Abstract. During the last decades, the global prevalence of dengue progressed 
dramatically. It is a disease which is now endemic in more than one hundred 
countries of Africa, America, Asia and the Western Pacific. This study ad¬ 
dresses a mathematical model for the dengue disease transmission and finding 
the most effective ways of controlling the disease. The model is described 
by a system of ordinary differential equations representing human and vector 
dynamics. Multiobjective optimization is applied to find the optimal control 
strategies, considering the simultaneous minimization of infected humans and 
costs due to insecticide application. The obtained results show that multiob¬ 
jective optimization is an effective tool for finding the optimal control. The set 
of trade-off solutions encompasses a whole range of optimal scenarios, provid¬ 
ing valuable information about the dynamics of infection transmissions. The 
results are discussed for different values of model parameters. 
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1. Introduction 

Dengue is a vector-borne disease transmitted from an infected human to a female 
Aedes mosquito by a bite. The mosquito, which needs regular meals of blood to 
feed their eggs, bites a potentially healthy human and transmits the disease, turn¬ 
ing it into a cycle. There are four distinct, but closely related, viruses that cause 
dengue. The four serotypes, named DEN-1 to DEN-4, belong to the Flavivirus 
family, but they are antigenically distinct. Recovery from infection by one serotype 
provides lifelong immunity against that serotype but provides only partial and tran¬ 
sient protection against subsequent infection by the other three viruses. There are 
strong evidences that a sequential infection increases the risk of developing dengue 
hemorrhagic fever. 

The spread of dengue is attributed to the geographic expansion of the mosquitoes 
responsible for the disease: Aedes aegypti and Aedes albopictus. The Aedes aegypti 
mosquito is a tropical and subtropical species widely distributed around the world, 
mostly between latitudes 35°N and 35°S. In urban areas, Aedes mosquitoes breed 
on water collections in artificial containers such as cans, plastic cups, used tires, 
broken bottles and flower pots. Due to its high interaction with humans and its 
urban behavior, the Aedes aegypti mosquito is considered the major responsible for 
the dengue transmission. The life cycle of a mosquito has four distinct stages: egg, 
larva, pupa and adult. In the case of Aedes aegypti, the first three stages take place 
in, or near, the water, whereas the air is the medium for the adult stage [Sl¬ 
it is very difficult to control or eliminate Aedes aegypti mosquitoes due to their 
resiliency, fast adaptation to changes in the environment and their ability to rapidly 
bounce back to initial numbers after disturbances resulting from natural phenomena 
(e.g., droughts) or human interventions (e.g., control measures). Primary preven¬ 
tion of dengue resides mainly in mosquito control. There are two primary methods: 
larval control and adult mosquito control, depending on the intended target. Lar- 
vicide treatment is done through long-lasting chemical in order to kill larvae and 
preferably have World Health Organization clearance for use in drinking water [3]. 
Adulticides is the most common measure. Its application can have a powerful im¬ 
pact on the abundance of adult mosquito vector. However, the efficacy is often 
constrained by the difficulty in achieving sufficiently high coverage of resting sur¬ 
faces [4]. 

The present study addresses a problem of finding the most effective ways of con¬ 
trolling the dengue disease. To this end, a mathematical model for the dengue 
transmission, including a control variable presented in |18) . is adopted. The main 
contributions are: (i) adapting multiobjective optimization methodologies for find¬ 
ing the optimal control in a mathematical model for the dengue transmission, (ii) 
analysis of the search ability of different methods on the resulting optimization 
problem, (iii) discussion of results for the optimal control problem, emphasizing 
advantages of the proposed approach when compared with results available in the 
literature. The aim of the work is to promote multiobjective optimization in epi¬ 
demiological research community and practice. 

The remainder of this paper is organized as follows. Section [2] describes a math¬ 
ematical model for the dengue disease transmission. The problem of finding the 
optimal control is formulated in Section [3l including traditional and proposed ap¬ 
proaches. Section 0] presents some general concepts in multiobjective optimization 
and methods adopted for solving the problem. Section [5] presents and discusses our 
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results, obtained by numerical simulations, which includes the comparison of differ¬ 
ent multiobjective optimization methods, control strategies and dengue dynamics 
that correspond to different optimal scenarios. The study is ended with Section |6] 
of conclusions and some directions of future work. 


2. The mathematical model 


This section introduces a mathematical model for the dengue disease transmis¬ 
sion. The model is based on a system of ordinary differential equations, and in¬ 
cludes the real data of a dengue disease outbreak that occurred in the Cape Verde 
archipelago in 2009 [HIIT]. 

The model consists of eight mutually-exclusive compartments representing the 
human and vector dynamics. It also includes a control parameter, an adulticide 
spray, as a measure to fight the disease. 

The notation used in the mathematical model includes four epidemiological states 
for humans: 

Sh{t) - number of susceptible individuals at time t; 

Eh(t) - number of exposed individuals at time t; 

Ih(t) — number of infected individuals at time t; 

Rh{t) — number of resistant individuals at time t. 

It is assumed that the total human population {Nh) is constant, so 

Nh = Shit) + Ehit) + Ihit) + Rhit). 


There are also four other state variables related to the mosquitoes (disease vectors): 

^m(f) “ number of vectors in aquatic phase at time t; 

Smit) - number of susceptible vectors at time t; 

Emit) - number of exposed vectors at time t] 

Imit) - number of infected vectors at time t. 

Similarly, it is assumed that the total adult mosquito population (AVi) is constant, 
which means that Nm = Smit) + Emit) + Imit)- 

The model includes a control variable, which represents the amount of insecticide 
that is continuously applied during a considered period, as a measure to fight the 
disease: 


c(t) - level of insecticide campaigns at time t. 


The control variable c(t) is an adimensional value that is considered in relative 
terms, varying from 0 to 1. 

In the following, for the sake of simplicity, the independent variable t is omitted 
when writing the dependent variables (for instance, Sh is used instead of Shit))- 
The parameters necessary to completely describe the model are presented in 
Table [TJ This set of parameters includes the real data related to the outbreak of 
dengue disease occurred in Cape Verde in 2009 m- 

Furthermore, in order to obtain a numerically stable problem, all the state vari¬ 
ables are normalized as follows: 


Sh = 


Nh' 


eh = 


Eh. 

Nh' 


ih = 


Nh' 


Th = 


Rh 

^h' 


Ah 

kNh 


Sr, 


mNh 


E„ 


mNh 


mNh 
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Parameter 

Description 

Value 

Nh 

total population 

480000 

B 

average daily bites (per mosquito per day) 

1 

^mh 

transmission probability from Im (per bite) 

0.375 

l^hm 

transmission probability from Ih (per bite) 

0.375 

d/ph 

average human lifespan (in days) 

1/(71 X 365) 

d/ih 

mean viremic period (in days) 

1/3 

1/ f-^m 

average lifespan of adult mosquitoes (in days) 

1/11 


number of eggs at each deposit per capita (per day) 

6 

PA 

natural mortality of larvae (per day) 

1/4 

VA 

rate of maturing from larvae to adult (per day) 

0.08 

1/ 

extrinsic incubation period (in days) 

1/11 

i/p/i 

intrinsic incubation period (in days) 

1/4 

m 

number of female mosquitoes per human 

6 

k 

number of larvae per human 

3 


Table 1. Model parameters. 


Thus, the dengue epidemic is modeled by the following nonlinear time-varying state 
equations [T8] : 

(^BPfjjihTtlyim “t“ ^h)Sh 

^ = BPmhmimSh - {vh + Hh)eh 
at 


( 1 ) 


dih 
dt 
drh 
dt 

- am)(Sm + em + im) “ {VA + 


= VhCh - ivh + ti'h)ih 

= Vhih - tJ-hrh 

m 

k 


dt 

dSm k ( TD R ' \ \ 

,, — VA \d> PfiTn^^h T pmj^m 

dt m 

dcui 

— dd (/4m “t“ CCm 

dim 
dt 


— ^m^m CXm 


subject to the initial conditions 

s^(0) = 0.99865, e,,(0) = 0.00035, 4(0) =0.001, rft(O) = 0, 
am(0) = 1, Sm(0) = 1, em(0)=0, i„(0)=0. 

Since any mathematical model is an abstraction of a complex natural system, 
additional assumptions are made to make the model mathematically treatable. This 
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is also the case for the above epidemiological model, which comprises the following 
assumptions: 

• the total human population (Nh) is constant; 

• there is no immigration of infected individuals into the human population; 

• the population is homogeneous, which means that every individual of a 
compartment is homogeneously mixed with other individuals; 

• the coefficient of transmission of the disease is fixed and does not vary 
seasonally; 

• both human and mosquitoes are assumed to be born susceptible, i.e., there 
is no natural protection; 

• there is no resistant phase for mosquitoes, due to their short lifetime. 


3. Problem formulation 


In this section, the problem of finding the optimal control in the mathematical 
model for the dengue transmission is formulated. The optimal control represents 
the most effective way of controlling the disease that can be adopted by author¬ 
ities in response to its outbreak. First, a traditional approach based on optimal 
control theory is presented. Then, the proposed approach, based on multiobjective 
optimization, is introduced. 

3.1. Optimal control theory. A problem of finding a control law for a given 
system is commonly formulated and solved using optimal control theory. Our aim 
is to find the optimal value c* of the control c, such that the associated state 
trajectories sj!^, e^, i'^, r^, a^, s^, e^, are solution of the system (HD in the 
time interval [0,T], subject to the initial conditions ([2]), and minimize an objective 
functional. 

Consider the state system of ordinary differential equations (HD and the set of 
admissible control functions given by: 

n = {c(-) G L“(0, T) I 0 < c(t) < 1, Vt G [0, T]}. 

The objective functional is defined by [l6] : 



( 3 ) 


where yu and ys are positive constants representing the costs weights of infected 
individuals and spraying campaigns, respectively. The objective functional given 
by dSD is a function of state and control variables, representing cumulative costs 
due to infected population and prevention measures. The optimal control problem 
consists in determining (s^, ejjj, r^, s^, e’^, i'^), associated to an admissible 

control c*(-) G on the time interval [0, T], satisfying ([T]), the initial conditions ([2]). 
and minimizing the cost functional ([3]), i.e.. 


J{c*{-))= min J(c(-)). 

c(.)6n 


The optimal control can be derived using Pontryagin’s maximum principle m- 

3.2. Multiobjective approach. An approach based on optimal control theory 
allows to obtain a single optimal solution. The obtained solution represents some 
decision maker’s perspective on controlling the disease. This is reflected by the 
constants yu and ys that must be provided in advance. Though the choice of 
proper values for the constants is not an easy task for all the cases. A single 
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optimal solution provides limited information for the decision maker, while leaving 
a large range of alternatives unexplored. To overcome these drawbacks, the present 
study decomposes the cost functional ([3]) into two separate objectives and seeks the 
optimal control, minimizing simultaneously the costs due to the infected population 
and the costs associated with insecticide. The resulting optimization problem is 
defined as follows: 


(4) 


minimize: 


/i(*/«(•)) = Jj ihit)dt 
/2(c(-)) = fo c{t)dt 


subject to: o and 


where T is a given period of time, /i and /2 represent the cost incurred in the 
form of infected population and the cost of applying insecticide for the period T, 
respectively. 


4. Multiobjective optimization 

This section presents general concepts in multiobjective optimization and dis¬ 
cusses some popular methods for multiobjective optimization based on scalariza- 
tion. The outline of an algorithm for finding optimal solutions is presented. For a 
method that gives a direct solution to the full Pareto set, via solutions of a certain 
PDF, we refer the reader to IS Eli. For other frameworks see m and references 
therein. 

4.1. General definitions. Without loss of generality, a multiobjective optimiza¬ 
tion problem (MOP) with m objectives and n decision variables can be formulated 
as follows: 

minimize: f (x) = {fi{x), f 2 ix),..., 

^ subject to: a; S 

where x is the decision vector, C R" is the feasible decision space, and f{x) is 
the objective vector defined in the attainable objective space 0 C R"*. 

In multiobjective optimization, the Pareto dominance relation is usually used to 
define the concepts of optimality. For two solutions x and y from fl, a solution 
X is said to dominate a solution y (denoted by a: ^ y) if for all z G {1, ..., m} 
Mx) < fi{y) and there exists at least one objective such that fj{x) < fj{y). 

A solution x* G il is Pareto optimal if and only if: 

$y G ■ y < X* . 

In the presence of multiple conflicting objectives, there is a set of optimal solu¬ 
tions, known as the Pareto optimal set. For MOP ((S]), the Pareto optimal set (or 
Pareto set for short) is defined as: 

rS = {x* Gn\$y Gil: y^ x*}. 

For MOP ([5]) and the Pareto set VS, the Pareto optimal front (or Pareto front 
for short) is defined as: 

VV={f{x*) G 0|a;* e VS}. 

Since it is often not possible to obtain the whole Pareto set, solving ([SJ is usually 
understood as approximating the Pareto set by obtaining a set of solutions that are 
as close as possible to the Pareto set and as diverse as possible. 
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4.2. Scalarization methods. In the following, four scalarization methods for mul¬ 
tiobjective optimization, able to deal with convex and nonconvex Pareto fronts, are 
discussed. 

4.2.1. The e-Constraint method. In the e-constraint method [5], one of the objective 
functions is selected to be minimized, whereas all the other functions are converted 
into constrains by setting an upper bound to each of them. It can be defined as: 

minimize: fi{x) 

subject to: fi{x) < e^, Vz € {1,..., m} A z ^ 

where the Ith objective is minimized and the parameter represents an upper 
bound of the value of fi. 

4.2.2. Chehyshev’s method. Chebyshev’s method [T] minimizes the weighted metric 
Lp, with p = oo, which measures the distance from any solution to some reference 
point. It can be formulated as: 

(6) minimize: max {wi{fi{x) — z*)}, 

l<2<m 

where z* = (z *,..., is a reference point and w = {wi ,..., Wm)"^ is a weight 

m 

vector such that Wj = 1. 

4.2.3. The goal attainment method. The problem shown in ([6]) can be reformulated 

as [TT] : 

minimize: a 

a3Gf2, Q;>0 

subject to: wi{fi{x) — z^) < a, 

^ Wmifmix) — < O', 

where z* = (z*,..., z^)"^ is a reference point, w = (wi,..., Wm)'^ is a weight vector 

m 

(J^wj = 1) and X G ft, a G K+ are variables. This method is often referred to as 
the goal attainment method m or the Pascoletti-Seraflni scalarization M- 

4.2.4. The normal constraint method. The normal constraint (NC) method [SJ [TU] 
minimizes one of the objective functions and uses an inequality constraint reduction 
of the feasible objective space. It can be formulated as: 

minimize: fi{x) 

subject to: vj{f{x) — z) <0, Vz G {1,...,m} Ai ^ I, 

where Vt is a vector from the zth corner, JT*, to the fixed ^th corner, Jif*, of the 
Pareto front 

Uj = 71 ^* — 71 **, Vz G {1,...,m} Ai ^ 1] 

2 : is a point on the hyperplane corresponding to a given vector w = {wi ,..., Wm)"^, 

z = ^w, 

m 

where $ = (7*^*, ■ • ■ ,T^*) is a m x m matrix and X) 
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To cope with differently scaled objectives, the NC method normalizes the objec¬ 
tive vector as: 


(7) 




fi - z'i 


ideal 


Vi G {1,... ,m}, 


where and are the ith components of the ideal and nadir points, respec¬ 
tively. 


4.3. An algorithm for the Pareto set approximation. The above discussed 
methods convert a multiobjective optimization problem into a single-objective prob¬ 
lem depending on some parameters. Solving the corresponding problem for different 
parameter settings allows to approximate multiple Pareto optimal solutions. The 
outline of an approach to approximate the Pareto set of problem (|3]) is shown in 
Algorithm [TJ 


Algorithm 1 Pareto set approximation 

1 : initialize: B 

2 : {}; 

3: for (3 G B do 

4: for initial point , find x* that minimizes fp{x,l3)] 

5: A-(-Aula:*}; 

6 : ^ X*] 

7: end for 
8: output: A; 


In line 1 of Algorithm[l] an initial point is generated as the null vector, x^^'> = 0 , 
and a set of parameters, B, used for scalarization, is initialized. For the e-constraint 
method, (3 G B corresponds to the value of e. In the case of the Chebyshev and 
goal attainment methods, /3 is a tuple of the form whereas /3 is a tuple 

of the form {ir, for the normal constraint method. In lines 3-7, for 

each corresponding scalarizing function, /^(a;,/3), a minimizer, a:*, is found using a 
single-objective optimizer and added to the approximation set, A. The minimizer of 
the previous scalarizing function becomes an initial point for solving the subsequent 
problem, to ensure efficiency of the approach. 

5. Numerical experiments 

This section provides results of the numerical experiments conducted for the 
dengue transmission model. The comparison of different scalarization methods is 
performed. The analysis of obtained solutions and the dengue dynamics are pre¬ 
sented. 

5.1. Experimental setup. The fourth-order Runge-Kutta method is used to dis¬ 
cretize the control and state variables of the system ([T]). The period [0,84] of 84 
days is discretized using the equally spaced time intervals of 0.25 (6 hours). This 
results in an optimization problem with x G [0,where x denotes the discrete 
control in O- The integrals defining the objective functionals in Q are calculated 
using the trapezoidal rule. 

For approximating the Pareto set of (j4]), scalarization is performed by defining 
and solving 100 scalarizing problems corresponding to the above discussed methods. 
As a single-objective optimizer (line 4 in Algorithm [IJ , the MATLAB® function 
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method 

measure 

e-Constraint 

Chebyshev 

Goal attainment 

Normal constraint 

Hyper volume 

0.972798 

0.959887 

0.959889 

0.977248 


Table 2. Hypervolume values for trade-off solutions obtained us¬ 
ing different scalarization methods (the higher the better). The 
best value is marked bold. 


fmincon with a sequential quadratic programming algorithm is used, setting the 
maximum number of function evaluations to 20, 000. 

5.2. Different scalarization methods. For comparison, the outcomes obtained 
by different scalarization methods are assessed using the hypervolume |21) . The hy¬ 
pervolume can be defined as the Lebesgue measure, A, of the union of hypercuboids 
in the objective space: 

HV = a( IJ [/i(a),ri] X • • ■ X [/„(a),r„] 

\aGA/\r 

where A = {ai,...,a|yi|} denotes a set of nondominated solutions and 

r = (ri,..., rm)"^ is a reference point. The hypervolume calculates a portion of the 
objective space dominated by A. It can measure both convergence and diversity. 
The higher the value of HV, the better the quality of A. For calculating the 
hypervolume, the objectives are normalized using 0 and [1,1] is used as a reference 
point. 

Table [5] presents the results for the four scalarization methods with respect to the 
hypervolume. The best result is obtained by the normal constraint method followed 
by the e-constraint method. Both methods minimize /i, when /2 is used to define 
the constraint for the objective space reduction. Since the NC method uses evenly 
distributed points on the ideal plane, the final points on the Pareto front approx¬ 
imation are likely to be more evenly distributed than using the usual e-constraint 
method. This allows to achieve the highest hypervolume. The Chebyshev method 
provides a slightly worse result than the goal attainment method. This can be be¬ 
cause the scalarizing function used by the Chebyshev method is nondifferentiable, 
which usually poses additional difficulties for optimization. 

Figure [T] shows the Pareto front approximations obtained by the methods. The 
plots confirm the previous observations based on the hypervolume values. The 
Chebyshev and goal attainment methods provide visually similar results. The e- 
constraint method attempts to find evenly distributed points with respect to the 
uniform division of / 2 , whereas the NC method seeks a uniform distribution of 
points, according to points on the hyperplane passing through the corner points of 
the Pareto front. Since the NC method gives clearly better results for o, solutions 
obtained by this method will be employed for the analysis of the dengue transmission 
in what follows. 

5.3. Dengue dynamics. Figures [Hand [3] present the trade-off solutions obtained 
for different values of transition probabilities Phm and Pmh, respectively. As it 
was seen during discussion of different scalarization methods, the trade-off curve 
for (jl]) exhibits a smooth relation between the insecticide cost, / 2 , and the infected 
human population, /i, in the range 0 < /2 < 5. This suggests that implementing 
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Figure 1. Trade-off curves obtained by different methods. 



(a) 0 < /2 < 10 (b) 10 < /2 < 90 


Figure 2. Trade-off curves for different values of Ptm- 


optimal control strategies in this range can produce the most signihcant reduction 
in infected individuals, having the highest ratio between the desirable effect and the 
cost. However, starting from some further point, reducing the number of infected 
humans can be achieved through exponential increase in spendings for insecticide. 
Even a small decrease corresponds to a high increase in expenses for insecticide. 
Scenarios represented by this part of the trade-off curve can be unacceptable from 
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(a) 0 < /2 < 10 (b) 10 < /2 < 90 


Figure 3. Trade-ofF curves for different values of pmh- 




(a) i;, for c(') = 0 and (b) for c(-) = 0 and 

varying /3/im varying 



(c) ih for c(-) = 1 and (d) for c(-) = 1 and 

varying d/im varying /3mh 


Figure 4. Dynamics of the state variable iu for varying values of 
Phm and Prnh, Corresponding to extreme solutions. 


the economical point of view. Furthermore, it should be noted that even with 
the maximum spending it is not possible to eradicate the disease, with the lowest 
obtained value of the infected population being 0.0042. 
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To provide a better visualization of the most interesting parts of the trade-off 
curves in Figures [2] and O the plots are separated into two parts and presented 
in different figures, according to the range of / 2 . From Figures [2a] and [3al it can 
be seen that the higher the corresponding transition probability, the larger number 
of infected individuals. Though this difference reduces as more control is applied, 
vanishing for the extreme scenario with c(-) = 1. Varying Phm and Pmh produces 
apparently similar effects on the shape of the Pareto front. From Figures [2a| and [5^ 
it can be observed that for higher values of Phm and I3mh the part of the trade-off 
curve in 0 < /2 < 5 becomes increasingly nonconvex. Solutions in this region can 
be unattainable for methods that face difficulties in dealing with nonconvexities in 
the Pareto front. 

Figure |4] illustrates dengue epidemics, corresponding to extreme scenarios (c(-) = 
0 and c(-) = 1). It can be seen that for higher values of transition probabilities the 
peak in infected population occurs early, having larger values (Figures |4a| and [4b|. 
On the other hand, applying the maximum control, c(-) = 1, the dynamics of ih{t) 
remain unchanged for different values of I3hm and f3mh (Figures 0^ and [4d|. 

To discuss intermediate scenarios, the concept of the knee solution of the Pareto 
front is adopted. According to [2], it can be defined as follows. Given a Pareto front 
with the normalized objectives, a boundary line T(_Pi,P 2 ) is constructed though 
two extreme points pi and P 2 - For any point z on the boundary line L{pi,p 2 ), a 
point pz on the Pareto front along the normal n of the boundary line is identified. 
The Pareto optimal point Pz» with the maximum distance from its corresponding 
boundary point z* along the normal direction is defined as the knee point. This 
concept is illustrated in Figure[6l The knee solution of the Pareto front is attractive, 
since it is often considered as the optima in objective trade-offs. 

Figure [5] shows the dynamics of the state and control variables for different values 
of the transition probabilities, corresponding to the knee solution. It can be seen 
that more controls are required for higher values of /3hm and Pmh during the entire 
period of study, which is consistent with expectations. However, a higher value of 
transition probability does not necessarily mean a higher value of infected human 
population for the knee solution. This is because the change in the shape of the 
Pareto front caused by varying a value of transition probability. 

From the Karush-Kuhn-Tucker condition, it can be induced that under smooth¬ 
ness assumptions the Pareto set of a continuous multiobjective optimization prob¬ 
lem defines a piecewise continuous (m — l)-dimensional manifold. The Pareto set 
of a continuous biobjective optimization problem is a piecewise continuous curve 
in R" [50]. The optimal control for (0|) can be defined as a surface in ih{t) x t. 
Solving (IH) gives a discrete representation of this surface. Thus, Figure |7| shows 
a discrete representation of a surface defining the Pareto set for ([4|). For a given 
value of /i, slicing the surface gives the optimal control trajectory. From the figure, 
it can bee seen how this dynamic changes from c(-) = 0 to c(-) = 1. A peak is 
observed in an early period of T, which grows in accordance with the decrease in 
the number of infected humans. Similarly, Figure |S| shows a discrete representation 
of a surface defining dynamics of infected humans across the Pareto set. A negative 
correlation between the optimal control and infected population can be observed. 
As the amount of control decreases, the number of infected humans increases. For 
higher values of the control, the peak in ih{t) is smaller and occurs later. When the 
control is reduced, the peak of disease outbreak grows up taking place earlier. 
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(a) c for knee point and (b) c for knee point and 

varying I3hm varying 0^^ 




t t 

(c) ih for knee point and (d) if^ for knee point and 

varying 0hm varying 0mh 


Figure 5. Dynamics of control and state variables for varying 
values of j3hm and /3mh, corresponding to knee solution. 



Figure 6. Knee point identification. 
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Figure 7. Discrete representation of the Pareto set. 




Figure 8. Discrete representation of infected humans. 


From the above discussion, it can be seen that the results for the optimal control 
problem, presented in Section [31 that are obtained using a multiobjective optimiza¬ 
tion approach provide comprehensive insights about optimal strategies for dealing 
with the dengue epidemic and dynamics resulting from implementing those strate¬ 
gies. The ability of trade-off solutions to reflect the underlying nature of the prob¬ 
lem constitutes a major advantage of multiobjective optimization, motivating its 
practical use in the process of planning intervention measures by health authorities. 

6. Conclusions 

Due to difficulties in treating the dengue disease, controlling and preventing 
its outbreaks is essential for keeping people healthy, especially in regions where 
the threat of dengue is high. This study discussed a mathematical model for the 
dengue disease transmission from the optimal control point of view. Multiobjective 
optimization approach is suggested for finding the optimal control strategies to deal 
with an outbreak of the dengue epidemic. A biobjective optimization problem is 
formulated, involving minimization of expenses due to the infected population and 
costs of applying insecticide. The approach avoids the use of a priori information 
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provided by the decision maker in the form of weight coefficients and allows to 
reflect the intrinsic nature of the problem. 

The problem is numerically solved by discretizing the control variable and using 
scalarization methods for approximating multiple Pareto optimal solutions. The 
obtained results reveal different perspectives on applying insecticide: a low number 
of infected humans can be achieved spending larger financial resources, whereas low 
spendings for prevention campaigns result in significant portions of the population 
affected by the disease. Different trade-offs between the objectives are represented 
by the obtained solutions. Varying transmission probabilities introduces changes 
in relation between the objectives and their trade-offs. Because the problem is 
convex, the Pareto set defining the optimal control remains unchanged for different 
values of transmission probabilities. The results of the study suggest advantages of 
multiobjective optimization for finding the optimal control. Once the Pareto set is 
approximated, the final decision on the control strategy can be made taking into 
consideration available financial resources and goals of public health care. 

As future work, it is intended to extend the model for the dengue transmission 
including different control variables that represent distinct measures for controlling 
the disease. Investigating impacts of different values of model parameters is also 
the subject of future work. 
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